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Times of arrival of high energy neutrinos encode information about their sources. We demonstrate 
that the energy-dependence of the onset time of neutrino emission in advancing relativistic jets 
can be used to extract important information about the supernova/gamma-ray burst progenitor 
structure. We examine this energy and time dependence for different supernova and gamma-ray 
burst progenitors, including red and blue supergiants, helium cores, Wolf-Rayet stars, and chemically 
homogeneous stars, with a variety of masses and metallicities. For choked jets, we calculate the 
cutoff of observable neutrino energies depending on the radius at which the jet is stalled. Further, 
we exhibit how such energy and time dependence may be used to identify and differentiate between 
progenitors, with as few as one or two observed events, under favorable conditions. 
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I. INTRODUCTION 

There is growing observational evidence and the- 
oretical foundation for the connection between core- 
collapse supernovae (CCSNe) and long gamma-ray bursts 
(GRBs) [1-3]. Location of GRBs in blue-luminosity re- 
gions of their host galaxies, where massive stars form 
and die, and CCSN signatures in the afterglows of nearby 
GRBs have provided strong evidence for the CCSN-GRB 
connection (see summaries by [4-6] ) . While many details 
of the CCSN-GRB relationship are still uncertain, cur- 
rent theoretical models suggest that a canonical GRB has 
a relativistic jet from a central engine [7-13]. 

Gamma rays are emitted by high energy electrons in 
the relativistic jet [14]. If the relativistic jet is able to es- 
cape from the star, the gamma rays can be observed and 
the GRB is coined "successful." If the jet stalls inside 
the star, however, the gamma-ray signal is unobservable 
and the GRB is coined "choked" [15]. While only - 1CT 3 
of CCSNe has extremely relativistic jets with V > 100, 
which lead to successful GRBs, a much larger subset of 
non-GRB CCSNe appears to be accompanied by colli- 
mated mildly relativistic jets with r ~ 10 [6, 16-20]. 
These CCSNe with mildly relativistic jets may make up 
a few percent of all CCSNe [17, 21-23]. Similarly, the 
jet's Lorentz factor for low-luminosity GRBs may also 
be much below those of high luminosity GRBs [24-27]. 

GRB progenitors may also produce high energy neu- 
trinos (HENs) [28]. The relativistic jet responsible for 
the gamma rays has been argued to be responsible for 
shock-acceleration of protons to ultrarelativistic energies, 
leading to nonthermal HENs produced in photomeson in- 
teractions of the accelerated protons. Recent upper lim- 
its from the IceCube detector [29] disfavor GRB fireball 
models with strong HEN emission associated with cos- 
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mic ray acceleration. However, milder HEN fluxes or 
alternative acceleration scenarios are not ruled out [30]. 
Moreover, the constraints weaken substantially when un- 
certainties in GRB astrophysics and inaccuracies in older 
calculations are taken into account, and the standard 
fireball picture remains viable [31, 32]. Proton accelera- 
tion may occur in external shocks [14, 33], both forward 
and reverse, as well as in internal shocks associated with 
the jet [34-40]. HENs, unlike gamma rays, can thus be 
emitted while the jet is still inside the star. Therefore, 
while gamma rays are expected to be emitted only if the 
GRB is successful, HENs are expected for both success- 
ful and choked GRBs. Thus, not only does the expected 
rate increase to a larger fraction of the CCSN rate of 
(2 - 3)yr~ 1 in the nearest lOMpc [41-43], these mildly 
relativistic jets also present a more baryon-rich environ- 
ment conducive to neutrino emission. Emission of HENs 
from internal shocks of mildly relativistic jets has there- 
fore been considered to be an important contribution to 
the overall observable neutrino flux [44-46]. 

As only neutrinos and gravitational waves can escape 
from the inner regions of a star, it means that they are 
a unique tool to study the internal structure of GRB 
progenitors. Whether and when HENs escape the star 
depend on the star's optical depth to neutrinos, which 
in turn depends on where neutrinos are emitted inside 
the star, as well as the neutrino energy - HENs can only 
escape after the relativistic jet reaches a radius where 
densities are low enough. This finite neutrino optical 
depth can delay and modify the spectrum of HEN emis- 
sion from GRBs. Detecting HENs with time and en- 
ergy information, at present and upcoming neutrino tele- 
scopes, therefore presents an unprecedented opportunity 
to probe GRB and CCSN progenitors, which may shed 
light on the GRB mechanism and the CCSN-GRB con- 
nection. 

Probing the interior structure of GRB/SN progeni- 
tors via HENs has first been suggested by Razzaque, 
Meszaros, and Waxman [47]. They showed, for two spe- 
cific stages of jet propagation, that the observable neu- 
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trino spectrum is affected by the stellar envelope above 
the jet head, which can in turn be used to examine this 
envelope via the detected neutrinos. Horiuchi and Ando 
[46] also mention HEN interactions for jets within the 
stellar envelope (see also Section IIIB). 

This article examines the question: What do the times 
of arrival of detected high energy neutrinos tell us about 
the properties of their source? We investigate the role 
of the opacity of CCSN/GRB progenitors in the proper- 
ties and distribution of observed HENs, and how these 
observed HEN properties can be used to probe the pro- 
genitors' structure. Studying the optical depth at which 
HENs can escape the progenitor, we find a progenitor- 
and energy-dependent temporal structure of the high en- 
ergy neutrino emission and jet breakout. Observations 
of HEN signatures of CCSNe or GRBs at neutrino tele- 
scopes, even with one or two events, could provide crucial 
information for differentiating between progenitors and 
characterizing their properties. Such information would 
advance our understanding of CCSNe, GRBs, and their 
relationship to each other. 

The paper is organized as follows: In Sec. II, we review 
CCSN and GRB models. In Sec. Ill, we briefly discuss 
HEN production in GRBs, propagation in the stellar ma- 
terial and flavor oscillations thereafter, and detection at 
Earth. In Sec. IV A, we describe our calculations of the 
neutrino interaction length and optical depth inside the 
stellar envelope, and present our results on the energy- 
dependent radius from which neutrinos can escape. In 
Sec. IV B, we discuss the temporal structure of energy- 
dependent neutrino emission from advancing and stalled 
jets for different stellar progenitors. This is followed in 
Sec. IV C by our results for energy dependent onset and 
emission duration of HENs. In Sec. V, we present our 
interpretations for the energy-dependent onset of HEN 
emission and discuss how it probes the progenitors, par- 
ticularly with a few detected neutrinos. We summarize 
our results in Sec. VI. 



II. CCSN AND GRB PROGENITORS 

Current understanding of canonical long GRBs sug- 
gests that they are collapsars requiring a massive pro- 
genitor star that is (i) rapidly spinning [48-51] and (ii) 
has a small radius (~ solar radius) [4, 52]. Successful 
GRBs also appear to prefer a lower metallicity [53], but 
choked GRBs may not require that. While this limited 
information does not always allow one to identify a spe- 
cific progenitor, it does suggest that the progenitors are 
massive rotating stars [4]. 

Rotating red and blue supergiants (e.g., [54, 55]) may 
be the progenitors of many GRBs. These stars are in the 
final stages of the pre-collapse evolution of massive stars, 
whose collapse can naturally lead to CCSNe and GRBs. 
Furthermore, some of these stars may lose their hydrogen 
envelope due to a binary companion, which can help the 
stars retain the fast rotation necessary for the creation 



of GRBs [56]. 

Wolf-Rayet (WR) stars are originally heavy, but lose 
their hydrogen envelope (and therefore a significant frac- 
tion of their mass) through stellar winds. A relativistic 
jet from a rotating WR star can therefore escape with- 
out having to penetrate a hydrogen envelope, making 
these stars a common type of progenitor [55, 57]. Mass 
loss through stellar winds is expected to be significant 
for stars with higher metallicity [58] . One difficulty with 
such mass loss is that it carries away crucial angular mo- 
mentum from the star. As the emergence of relativistic 
jets requires a very rapidly rotating core, losing angular 
momentum decreases the possibility of a GRB [59] . As a 
result, compact progenitors (i.e. that have lost their hy- 
drogen and/or helium envelopes) also regularly explode 
as type Ibc supernovae without indications of a central 
engine injecting jet power into the explosion (e.g., [60]). 

Alternatively to stellar winds, massive stars can lose 
their hydrogen envelope to a companion star, which can 
leave more angular momentum at the core. Such rotating 
objects, composed of the bare helium core left behind, 
may also be GRB progenitors. Unusually rapid rotation 
on the main sequence can also result in mass loss [59]. 

Single stars with extremely rapid rotation may experi- 
ence almost complete mixing on the main sequence [59] , 
leading to a chemically homogeneous star. Such stars 
bypass the red giant phase and resemble WR stars, but 
with little mass loss. This scenario is particularly inter- 
esting for GRB production as it combines low mass loss 
with rapid core rotation, the two prerequisites for GRB 
emission. 

Guided by these facts, we shall examine the following 
progenitor models in detail: 

1. Red supergiant - zero-age main sequence (ZAMS) 
mass of 15 M Q , both with solar (15S) and low (15L) 
metallicities. We further study the 15L model with 
its hydrogen envelope removed (15Lc) due to, e.g., 
a companion star. 

2. Wolf-Rayet star - ZAMS mass of 75 M with solar 
metallicity (75S). 

3. Bare helium core - ZAMS mass of 40 M with solar 
metallicity (40S). 

4. Chemically homogeneous star - ZAMS mass of 
16 M with low metallicity (16T). 

We have indicated, in parentheses, the names of the mod- 
els we consider in this study. To obtain the matter dis- 
tribution and composition of these models, we use the 
numerical results of Woosley et al. [55], Woosley and 
Heger [59], and Heger et al. [61]. A detailed list of their 
properties (and a reference to the literature) is given in 
the upper box in Table I. We also indicate in the ta- 
ble the references to the numerical results of the stellar 
progenitor models. 
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Model 


Z 


M 


M He 


M H 


R 




Ref. 




[Zq] 


[M ] 


[M s ] 


[M ] 


[10 13 cm] 


[10 11 cm] 




15L 


10" 4 


14.9 


4.5 


7.4 


0.4 


1.4 


[55] 


15Lc 


1(T 4 


5.2 


1.7 


0.4 


N/A 


1.4 


[55] 


15S 


1 


12.6 


4.0 


5.5 


5.9 


1.9 


[55] 


16T 


10~ 2 


15.1 


0.1 


3 x 10" 4 


8 x 10" 2 


N/A 


[59] 


40S 


1 


8.7 


0.1 


10" 3 


8 x 10" 3 


0.8 


[55] 


75S 


1 


6.3 


0.2 


8 x 10~ 4 


7 x 10" 3 


0.7 


[55] 


12L 


icr 4 


11.9 


3.4 


6.2 


2.4 


1.6 


[55] 


20L 


10~ 4 


19.9 


6.1 


9.0 


0.3 


1.9 


[55] 


25L 


10~ 4 


24.9 


7.7 


10.4 


0.3 


2.3 


[55] 


35L 


10~ 4 


34.8 


10.9 


12.7 


1.2 


3.4 


[55] 


12S 


i 


10.7 


3.4 


5.1 


4.3 


1.9 


[61] 


15Sb 


i 


11.9 


3.8 


4.8 


6.1 


2.3 


[61] 


20S 


i 


12.7 


4.0 


3.5 


7.7 


2.8 


[61] 


25S 


i 


12.2 


3.2 


1.7 


8.2 


3.2 


[61] 


35S 


i 


14.6 


2.2 


0.4 


0.1 


0.2 


[61] 



TABLE I. Properties of pre-supernova stellar models used in 
the analysis. The columns are: model name, metallicity, pre- 
supernova (PS) stellar mass, helium mass, hydrogen mass, 
stellar radius, helium core radius, and reference to models. 
Model names contain the ZAMS stellar mass, and a letter 
representing metallicity (L - low, T - 1% Solar, and S - Solar). 
There is additional differentiation between the various 15 M 
models we use. 



In addition to the models listed above, we examine 
other low- and solar-metallicity stars in the ZAMS mass 
range of (12 - 35) M Q , also listed in the lower box in 
Table I, to investigate some of our results' dependence on 
stellar mass and metallicity. Some of these models, which 
have ZAMS masses of < 30 M will probably not create 
successful GRBs. Nevertheless, they can have choked rel- 
ativistic jet activity that can be observed through neu- 
trinos. 



III. HIGH-ENERGY NEUTRINOS 
A. Production 

The variability in the output of the GRB's central en- 
gine results in internal shocks within the jet, which ac- 
celerate electrons and protons to high energies. Internal 
shocks can occur even when the relativistic jet is still 
propagating inside the star [15]. For jets inside the star, 
reverse shocks can also occur [46] at the head of the jet, 
which can also accelerate electrons and protons to rela- 
tivistic energies. 

Relativistic electrons emit gamma rays through syn- 
chrotron or inverse-Compton radiation. Relativistic pro- 
tons interact with these gamma rays (pry), or with other 
non- relativistic protons (pp), producing pions and kaons. 
Photomeson interactions produce charged pions (7r + ) 
through the leptonic decay p7 — > tt + . Proton-proton in- 



teractions produce charged pions (ir^) and kaons (if ). 
Charged pions and kaons from these processes decay into 
neutrinos through 

TT^i^ ->M ± +^(17 AI ) (1) 

Muons further decay to produce secondary neutrinos 
through, e.g., the process /i + — > e + +i/ e +V fJ- . However, if 
the synchrotron photon density is high enough, or in the 
presence of strong magnetic fields (i.e. for smaller radii), 
they may immediately undergo radiative cooling, giving 
a flux of lower energy neutrinos [44, 45]. 

The energy of charged mesons from both p7 and pp 
processes is about 20% of the proton's energy, while 
roughly 1/4" 1 of this energy is given to v^ip^) [36]. The 
energy of the produced neutrinos and antineutrinos is 
therefore ~ 5% of the proton energy. The energies of 
the photon (e 7 ) and proton (e p ) in the pry interaction 
need to satisfy the photo-meson threshold condition of 
the A-resonance [36] 

e 7 e p w 0.2r 2 GeV 2 (2) 

where T is the Lorentz factor of the shock. Assuming 
Lorentz factor T ~ 300 and observed 7 ray energy e 7 ~ 
1 MeV, one obtains a characteristic neutrino energy of 
e v ~ 10 14 eV [36]. 

For both pj and pp processes, the energy spectrum of 
HENs is determined mainly by the proton energy spec- 
trum, and the optical depths of the pry and pp interac- 
tions. The distribution of the proton energy e p in inter- 
nal shocks, in the observer frame, is d 2 N/(de p dt) cx e~ 2 , 
with a maximum energy cutoff due to photo-pion losses. 
The cutoff energy depends on both jet properties and the 
radius where the internal shocks occur. 

The collisional or radiative nature of the internal 
shocks at high densities can be an impediment to Fermi 
acceleration. For mildly relativistic jets, particle acceler- 
ation happens less efficiently because the shock is some- 
what spread out [62]. However, the details of the acceler- 
ation process are still uncertain, and for many alternative 
acceleration scenarios this is not an issue [63-65] . On the 
other hand, neutrinos are the best probes for the physics 
responsible for the acceleration, and HEN observation 
should shed light on this aspect. 

B. Interactions and Oscillations 

HEN emission from internal shocks can commence at 
a distance r s sw r 2 c<5i sw 3 x 10 9 cm [15], where Tj 
1 () is the jet Lorentz factor and St = 10~ 3 s is the jet 
variability time. Photo-pion losses determine the cut-off 
in the energy spectrum. 

Whether and when HENs escape the star depends on 
the star's optical depth to neutrinos, which depends on 
neutrino energy as well as where neutrinos are emitted 
inside the star. HENs can only escape after the relativis- 
tic jet reaches a low density region, beyond which the 
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interaction of HENs with the stellar medium is negligi- 
ble. 

The effect of HEN interaction with matter on observ- 
able HEN emission prior to the outbreak of the jet has 
been discussed previously, e.g., by Razzaque, Meszaros, 
and Waxman [47]. Using a simplified model for the jet 
and progenitor star, Razzaque et al. found that HEN 
interaction is negligible if the jet gets close to the surface 
(estimated as having roughly ~ 0.1 M© envelope mass 
over 47T above the jet front). If the jet is deeper inside 
the star (with overlying envelope material of ~ 1M©), 
they found HEN interaction effects to be noticeable, es- 
pecially for stellar models that lost their hydrogen en- 
velope. For this latter case, they found that the neu- 
trino optical depth becomes larger than unity for neu- 
trino energies e„ > 2.5 x 10 5 GeV. Horiuchi and Ando 
[46] also mention HEN interactions for jets within the 
stellar envelope. They find that only neutrinos with en- 
ergies less than e v < 10 2 GeV can escape a progenitor 
star from r ss 10 10 cm. 

The interaction probability of HENs increases with 
their energy, and can be non- negligible if a neutrino beam 
were to travel through large quantities of dense matter, 
as is the case for neutrinos produced inside massive stars. 
The mean free path of HENs is determined by inelastic 
scattering processes, with neutrino-nucleon interactions 
(uN or VN) being the major determinant of the opti- 
cal depth [66, 67]. At the relevant neutrino energies, 
the interaction cross sections are approximately the same 
for all neutrino flavors [67], therefore we treat the cross 
section to be flavor independent (see [67]). Electron- 
antineutrinos are the exception to the above, as the inter- 
action of electron-antineutrinos and electrons (V e e) be- 
comes the dominant effect around the resonant neutrino 
energy e Ve 6.3 x 10 6 GeV [66]. We shall neglect this 
exception as it only results in the attenuation of a small 
fraction of HENs, i.e. the radius at which neutrinos of a 
given energy can first escape does not change. 

To obtain HEN interaction lengths, we adopt the neu- 
trino cross sections obtained by Gandhi et al. [66] (other 
calculations give similar results; see, e.g., [68]). We ad- 
ditionally take into account nuclear effects (i.e. that 
both free and bound nucleons are present) calculated 
by Pena et al. [69]. Gandhi et al. [66] calculated the 
cross sections for charged and neutral currents, the sum 
of these two giving the total cross section a v . For the 
range of neutrino energies of interest here, neutrino cross 
sections are, to a good approximation, flavor invariant 
[66, 67]. We note that the cross sections for neutrinos 
and antincutrinos are somewhat different. 

The neutrino-nucleon interaction cross section <7„(e„), 
obtained numerically by Gandhi et al. [66], assumes 
that matter is in the form of free nucleons. To take 
into account the presence of bound nucleons (mostly 
helium), we approximate nuclear effects following Cas- 
tro Pena et al. [69]. The ratio of the total cross section 
cr v (A) corrected for nuclear effects over the cross section 
a v without correction (free-nucleon case) decreases with 



energy and with atomic mass number A. The effect is 
practically negligible (order of a few percent) below a 
neutrino energy of e v ~ 10 5 GeV. Above e„ ss 10 5 GeV, 
we approximate the energy dependence presented in [69] 
(sec Figure 2 therein) with the empirical function 

a„(e v ) ~Vl0 5 GcV/ ' [ > 

Given the neutrino interaction cross sections, one can 
obtain the interaction length (or mean free path) \ v via 

I r \ = ^p{r)w A {r)N m a v {e v ,A) , (4) 

^v\^Ui") a 

where p(r) is the stellar density as a function of the radial 
distance r from the center of the star, w^(r) is the mass 
fraction of elements with mass number A, N av is the Avo- 
gadro constant, and a v {e v ) is the vN interaction cross 
section as a function of neutrino energy e„ and atomic 
mass number A. For antineutrino interaction length Aj7, 
the difference compared to \ u is the energy-dependent 
ratio of neutrino and antineutrino cross sections. 

Besides scattering, neutrinos also undergo flavor mix- 
ing due to neutrino flavor oscillations. After the neu- 
trinos escape the star, the neutrinos travel a very long 
distance in space wherein the wave-packets of each mass 
eigenstate must separate. Therefore, neutrinos observed 
at Earth should be considered as incoherent superposi- 
tions of these mass-cigenstatcs. If the expected neutrino 
flavor ratio leaving the source is cj) Ve :<j) v :(j> Ur = 1:2:0, neu- 
trino oscillations from the source to the detector trans- 
form these ratios to 1:1:1 [70]. If muon radiative cool- 
ing is significant, then only muon neutrinos are produced 
through the decay of kaons and pions, i.e. the flavor ratio 
at the source will be 0:1:0. 

These ratios at the source and at the detector are 0:1:0 
and 1:2:2, respectively [70]. Other flavor ratios are possi- 
ble and have been explored [71]. As the energy and time 
profile of neutrino emission remains almost unchanged by 
oscillations, we do not consider these oscillation effects. 
We note that these ratios are modified by neutrinos oscil- 
lations inside the star that may be important for neutrino 
energies < 10 4 GeV [67]. 



C. Detection 

HENs traveling through Earth interact with the sur- 
rounding matter creating secondary particles, mostly 
muons. Cherenkov radiation from these muons is de- 
tected by neutrino detectors, which reconstruct direction 
and energy based on the detected photons [72]. While 
the IceCube detector [29] has a threshold of ~ 10 2 GeV, 
Earth starts to become opaque to neutrinos with energies 
e v ^ 10 6 GeV, decreasing the neutrino flux that reaches 
the vicinity of the detector after crossing Earth [73] . So- 
called upgoing neutrinos (i.e. neutrinos whose trajectory 
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crosses Earth before reaching the detector) with ener- 
gies above ~ 10 6 GeV are practically undetectable, as 
most of them are absorbed before reaching the detector 
(note that this is true for muon-neutrinos; tau-neutrinos 
can penetrate Earth even at higher energies [73]). For 
upgoing muons this sets the observable energy window 
at (10 2 — 10 6 ) GeV. Downgoing and horizontal neutri- 
nos (i.e. neutrinos whose trajectory reaches the detector 
without crossing Earth) are detectable at practically any 
energy [73]. The disadvantage of such directions is the 
much higher background noise from atmospheric muons 
that are, for these directions, not filtered out by Earth. 
Due to this high background, most HEN analyses only 
consider upgoing HEN events. The only exception is 
extremely- high energy neutrinos (> 10 6 GeV), as atmo- 
spheric muons seldom reach such high energies. 

There are several currently operating HEN telescopes, 
e.g., IceCube [29], a km 3 detector at the South Pole, and 
Antares [74] in the Mediterranean sea. Both Antares 
and a third detector at the lake Baikal [75] are planned 
to be upgraded into km 3 telescopes [76] . These km 3 -scale 
detectors observe HEN events and measure their energy, 
time, and arrival directions. While uncertainties remain 
in the emission mechanism and expected source flux [30] , 
depending on the model, these km 3 -scale detectors may 
observe (1 — 100) HEN events in the energy range (10 2 — 
10 6 ) GeV, for a typical GRB at 10 Mpc. 

The uncertainty of muon energy reconstruction in 
neutrino telescopes is ~ 0.3 in log 10 (e„/GeV) [77-79], 
becoming somewhat larger [79] for partially contained 
muons of greater energies. Timing resolution is expected 
to be on the nanosecond level [80, 81]. Angular pointing 
is expected to be < 1° at these energies, which may allow 
identification of the source independently. Although we 
will not perform a detailed simulation of the observable 
signal, we shall keep in mind these experimental param- 
eters. 



IV. CALCULATIONS AND RESULTS 

In this section, we characterize the effect of neutrino 
interactions on the observable HEN energies from jets. 
This effect can modify the neutrino flux from choked jets 
as well as from precursor neutrino emission for successful 
jets. We consider internal shocks for the calculation of 
the temporal structure and derive the energy-dependent 
onset time and duration for HEN emission. 



A. Optical Depth for High-Energy Neutrinos 

Given the mass distribution p(r) in a star, one can 
employ the expression for the neutrino mean free path 
described in the previous section to calculate the HEN 
optical depth of the star for a given distance from the cen- 
ter. We are interested in the innermost radius at which 
neutrinos can escape from the star. 
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FIG. 1. (top) High energy neutrino optical depths as a func- 
tion of distance from the center of the star for different neu- 
trino energies for a ZAMS 15 M© star with low metallicity. A 
vertical dashed line shows the radius of the helium core. As 
a comparison, we show the optical depth of the helium-core- 
only case indicated with thick lines. The horizontal dashed 
line shows r = 1. Above this line, the stellar envelope is 
opaque to neutrinos, (bottom) The star's mass density as a 
function of distance from the center. 



The optical depth r of the star at a distance ro from 
its center, towards neutrinos that are produced at ro and 
are moving radially outward, is 



Tw(e„,ro) 



1 



(5) 



where R is the stellar radius. We approximate neu- 
trino absorption such that neutrinos with energy e v can- 
not escape from below a critical radius r v for which 
T v (e u ,r u ) f» 1. 

We calculated the optical depths for the considered 
massive stellar models (see Table I) as functions of ra- 
dial distance r from the center of the star, as well as 
neutrino energy t v . In Figure 1, we exhibit the represen- 
tative behavior of the neutrino optical depth as a func- 
tion of distance from the center for a stellar model with 
M = 15 M Q ZAMS mass and low metallicity. We can see 
that the hydrogen envelope of the ZAMS M = 15 M 
star with low metallicity, for most relevant energies, is 
transparent to neutrinos. We also see that, for this stel- 
lar model, the helium core becomes opaque to neutrinos 
with all depicted energies around ss 10 10 cm. 

To obtain Ty CIS Bi function of neutrino energy, we in- 
verted Equation 5 to derive the critical radii. In Figure 2, 
we present the critical radii for neutrinos. In the lower 
panel, we present the ratio of critical radius of antineu- 
trinos and neutrinos (for the representative M = 20 M 
low- metallicity case). One can see that the neutrino and 
antineutrinos have very similar critical radii. The maxi- 
mum difference between the two radii is about 25%, in- 
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FIG. 2. (top) Critical radius r„ as a function of neutrino 
energy for the considered progenitor models (see Table I). The 
horizontal dashed line shows the pre-supernova helium-core 
radius for the ZAMS 15 M© star with low metallicity (model 
15Lc). (bottom) Ratio of cross sections for neutrinos and 
antineutrinos o„/<j v (data taken from [66]), and the obtained 
ratio of the critical radii for neutrinos and antineutrinos, as 
the function of neutrino energy, for the ZAMS 20 Mq star 
with low metallicity. 



dicating that our results on critical radii are valid for 
antineutrinos as well. 

We have investigated the dependence of the above re- 
sults on the mass and metallicity of the progenitor. To 
do so, we repeated the above exercise for models with low 
and solar metallicity, from (12-35) M listed in Table I, 
and plotted the critical radius as a function of energy 
and ZAMS mass, separately for low and solar metallic- 
ity. A contour plot indicating the effect of ZAMS stel- 
lar mass on r v for different metallicities and neutrino 
energies is shown in Figure 3. The critical radius that 
corresponds to the supernova-progenitor's helium core is 
indicated with dashed horizontal lines (if all neutrinos 
with e„ < 10 6 GeV can escape the core, the horizontal 
dashed line lies above the shown parameter space; this is 
indicated with arrows pointing upwards). 

Figure 3 shows that the critical radius lies within the 
stellar helium core for neutrino energies of e v < 10 5 GeV 
for practically all massive progenitors. Consequently, the 
most relevant neutrino energy range for observations will 
become observable somewhere beneath the stellar core, 
making neutrino observations from choked-GRBs rele- 
vant. The results further indicate that r v (e v ) is deter- 
mined predominantly by the pre-supernova mass (and 
not the ZAMS mass) of the progenitors. 

Beyond the optical depth of the stellar envelope, here 
we examine the possibility of neutrino absorption by the 
jet itself and the jet head. For the fiducial values of 
our model and the shock radius (see below; luminosity 
Li SO = 10 52 ergs -1 , jet gamma factor Fj — 10, and shock 



radius r s — 3 x 10 9 cm), the density of the jet [see Eq. 
(8)] at r s in the observer frame is p w 0.3 gem -3 . This 
density is negligible compared to the stellar density p(r = 
3 x 10 9 cm) «3x 10 3 gcm -3 (model 15L), therefore the 
jet itself will not play a role in neutrino absorption. 

To estimate the density of the jet head, we take stellar 
density p(r = 3 x 10 9 cm) «3x 10 3 gcm~ 3 (model 15L), 



1/4-1/2 1/4 
52 ' 10 



and jet head Lorentz factor Fh = 0.1L 
0.13 for the fiducial values [15]. While the density of the 
jet head is significantly smaller than the stellar density, 
the shocked stellar density will be a few times greater 
than the stellar density [46, 82]. Nevertheless, the thick- 
ness of this shocked region is small compared to the stel- 
lar radius. Taking the simulation of Aloy et al. [82] 
as an example, the shocked stellar region at its peak is 
~ 2.5 x greater than the stellar density when the jet head 
is at 6 x 10 9 cm, and the thickness of the jet head is 
8 x 10 8 cm. With such jet head at 6 x 10 9 cm, the stel- 
lar column depth increases by ~ 20% (using model 15L). 
While such shocked stellar region does not change neu- 
trino absorption substantially compared to the rest of the 
star, it further increases the role of neutrino absorption 
in the observable neutrino spectrum. Further, for faster 
advancing jets, the jet carries more matter in front of it 
as it takes time for matter to flow sideways [83], making 
the shocked head even denser (see also [46]). For simplic- 
ity, in the results below we neglect the absorption due to 
the shock region. 



B. Temporal Structure of Jets 

Jets have been studied numerically and analytically, 
and a detailed understanding has been developed; see, 
e.g., references in [46]. For our work, we use the semi- 
analytic method of Horiuchi and Ando [46] to calculate 
the velocity of the jet head advancing inside the star 
in order to characterize the temporal structure of HEN 
emission. This simple treatment assumes a constant jet 
opening angle, and will suffice to illustrate our point, but 
in the future this may be improved through more detailed 
numerical modeling of the jet morphology for individual 
progenitors. 

Horiuchi and Ando consider the propagation of a rel- 
ativistic jet with Lorentz factor Tj ^S> 1. As the head of 
the jet advances through stellar matter with Lorentz fac- 
tor r^, a reverse and forward shock occur. The reverse 
shock decelerates the head, while a forward shock accel- 
erates the stellar material to Fh- In the following, we use 
the subscripts j (jet, unshocked), h (jet head, shocked), s 
(stellar, shocked), and ext (stellar, unshocked) to denote 
quantities at different regions in and around the jet. Us- 
ing this notation, the evolutions of the two shocks at the 
jet head are governed by the following equations [84, 85]: 

e s /n s m p c 2 = F h - 1, n s /n ext = 4F h + 3, (6) 
e h /n h m p c 2 = F h - 1, n h /n 3 = 4F h + 3, (7) 
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FIG. 3. Critical radius r v — ro(r v = 1) as the function of ZAMS stellar mass and neutrino energy e„ for low-metallicity (left) 
and Solar- met allicity (right) simulations by Woosley et al. [55], and Heger et al. [61]. The horizontal dashed lines show, for 
each ZAMS mass, the energy for which the critical radius is equal to the pre-supernova helium core radius (if all neutrinos 
with e v < 10 6 GeV can escape the core, the horizontal dashed line lies above the shown parameter space; this is indicated with 
arrows pointing upwards). Consequently, neutrinos with energies above the dashed lines cannot escape if produced inside the 
helium core. Note that, as shown for an example in Figure 1, the presence of the hydrogen envelope has only small effect on 
this threshold energy. 



where m p is the proton mass, c is the speed of light, and 
rii and e$ are the particle density and internal energy 
measured in the fluids' rest frames. Lorentz factors are 
measured in the lab frame, except Th = TjTh(l — /3jj3h) 
that is measured in the jet's comoving frame (fiiC is the 
velocity). For a jet with constant opening angle, the jet 
particle density at radius r in the jet frame is 



n 3 (r) 



4irr 2 T'jmpC 3 



(8) 



where Lj so is the isotropic-equivalent jet luminosity. At 
the jet head, the shocked jet head and shocked stellar 
matter are separated by a contact discontinuity, and are 
in pressure balance. Equating their pressures (ph = e/j/3 
and p s — e s /3, respectively) and using Equations (7) 
and (8), we arrive at 

n 3 _ (4T h + 3)(r h -l) 



(9) 



(4T h + 3)(r h -l) 

We numerically solve Equation (9) using the appropri- 
ate stellar particle number density n ext = n ext (r) of the 
considered stellar models (see Table I) to obtain the jet 
propagation velocity as a function of radius. We took 
the terminal Lorentz factor of the jet head to be Tj in 
the limit of zero density. Similarly to Ref. [46], we find 
that the velocity of the jet head is practically indepen- 
dent of the gamma factor Tj of the jet, and increases 
with the isotropic equivalent energy Li SO of the jet (in 
the relativistic limit, the gamma factor of the jet head 



r oc r 



1/4 



[15]). 



C. Energy-dependent Emission Onset Time 

The skewed neutrino emission due to stellar neutrino- 
opacity can be characterized through the temporal struc- 
ture of observable emission. Let to be the time when 
the relativistic jet reaches the shock radius r s sa r 2 c<5i « 
3x 10 9 cm. At this radius, internal shocks can form, caus- 
ing HEN production to commence [36]. Neutrinos that 
are produced while the jet is still beneath the stellar sur- 
face may only be able to escape from the star through 
the envelope at a later, energy- dependent "escape" time 
t e (e v ), where e„ is the energy of the neutrino. At t e (e v ) 
the jet has advanced far enough so the remaining enve- 
lope and the jet itself is no longer opaque to neutrinos 
with energy e„. More specifically, we take t e (e v ) to be 
the time when the jet reached the critical radius r v for 
which r„(e„, r„) = 1. 

In Figure 4 (left), we show t e (e v ) — to (observer frame) 
as a function of e„ for different stellar progenitors, which 
we calculated assuming a mildly relativistic (Tj = 10) 
jet with Li SO — 10 52 ergs -1 output. We can see that as 
neutrinos with lower energies can easily escape through 
the stellar envelope, one finds t e (e v <C 100 GeV) — to ~ 0. 

It is also interesting to compare the escape time t e (e v ) 
with the jet breakout time tb r from the stellar envelope. 
The point of jet breakout was chosen to coincide with the 
jet reaching the radius at which the simulated stellar pro- 
genitor models end. This corresponds to a sharp drop in 
matter density, dropping below 10 -10 gcm -3 for models 
15L and 16T, and dropping below 10 _4 gcm -3 for mod- 




FIG. 4. (left) Onset of observable HEN emission measured from the time when HEN production commences (t e (e v ) — to), as a 
function of neutrino energy, for different stellar progenitors, (right) Time of jet breakout measured from the time of the onset 
of observable HEN emission (tbr — t e (e v )), as a function of neutrino energy, for different stellar progenitors. The calculations 
are carried out for the stellar models in Table I with jet energy Li ao — 10 52 ergs -1 and jet Lorentz factor Tj = 10. 



els 40S and 75S, while the simulation ends at 10 13 cm for 
model 15S. Due to the low densities at the boundary of 
the simulated progenitors, the results should be robust to 
the specific choice of jet breakout radius. We note that 
the bulk of the gamma-ray emission from the jet may be 
shortly delayed compared to the breakout due to dissipa- 
tion within the jet until the jet advances to a distance of 
> 10 13 cm [13, 86-88]. The GRB will become detectable 
within a few seconds after the jet head leaves the helium 
core [89]. 



Figure 4 (right) shows the time it takes for the jet to 
break out of the star from the point from which neutri- 
nos can first escape (i.e. tbr ~ t e {^v)), as a function of e v 
for different stellar progenitors. As before, a mildly rela- 
tivistic (Tj = 10) jet with Li SO = 10 52 ergs -1 output has 
been used. A disadvantage of using tb r for comparison is 
that it is only available for successful jets, and therefore 
it cannot be used to characterize, e.g., choked GRBs. 



As the figure shows, when the jet is close to the center 
of the progenitor, only low energies neutrinos can escape 
the star. Gradually, as the jet proceeds outward, higher 
energy neutrinos become observable. Consequently, neu- 
trinos escaping (and observed) earlier are expected to 
have lower energies, leading to a time dependent neutrino 
energy distribution, with the average energy increasing 
with time. Such distribution may be indicated upon 
the detection of multiple neutrinos from a CCSN. The 
precise relationship of the energy-dependent onset times 
and emission durations of HENs can encode information 
about the progenitors' density profile and composition. 



D. Dependence on source parameters 

The results presented above were calculated for sources 
with mildly relativistic jets (Tj = 10) with Li SO luminos- 
ity, and with jet variability of St. From these parameters, 
Tj and St determine the shock radius r s (see Sec. IV C). 
For greater Tj and St, neutrino production will commence 
only at a greater radius, changing the lowest radius and 
earliest time from which neutrinos are observable. The 
jet luminosity Li SO affects the jet head velocity, with 
greater (smaller) Lj so corresponding to greater (smaller) 
velocity, which in turn change the time scale on Figure 4; 
the jet head velocity is practically independent of Tj and 
St (see also [46]). These dependencies, however, do not 
affect our conclusions qualitatively. 



V. INTERPRETATION OF ENERGY AND 
TIME STRUCTURE OF NEUTRINO EMISSION 

A. Strong-signal limit 

For a very nearby CCSN or GRB, one expects a mod- 
estly large HEN flux of perhaps > 100 events. Although 
this possibility is rare, it does represent a rather lucra- 
tive opportunity. The highest energy neutrino detected 
from such a source at any time constrains the total mat- 
ter content of the envelope above the jet head at that 
time. Therefore, one can take the opportunity to map 
the inner stellar density profile and the jet's velocity. 

Comparison of the emission onset profiles of different 
stellar progenitor models in Figure 4 (left panel) shows 
that, with a high-enough HEN flux, the stellar models 
have distinguishable energy-time profiles (i.e. cutoff ener- 
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gies as a function of time). For example, for a 5s long pre- 
cursor emission, the groups of models {16T}, {40S,75S} 
and {15L, 15Lc,15S} are distinguishable, while members 
of a given group have practically the same energy cutoff. 

Another consequence of the radius-dependent neutrino 
energy cut-off is that neutrinos emitted by choked rela- 
tivistic jets (e.g., choked GRBs) will have an energy cut- 
off, resulting in lower average energy than jets that suc- 
cessfully break through to the surface of the star (e.g., 
successful GRBs). Such a difference can possibly be in- 
dicated upon the detection of a sufficient number of neu- 
trinos from a set of CCSN with and without electromag- 
netic counterparts. The indication of a difference in av- 
erage energies from successful and choked jets can pro- 
vide information, e.g., on the distance the choked jets 
advance before they stall, a possible indicator of the en- 
ergy and/or baryon content of the jets. 



B. Weak-signal limit 

Most observed astrophysical sources will only lead to 
0(1) detected neutrinos. Therefore, we find it more prac- 
tical to consider the possibility that only a few neutrinos 
are detected. In the following, we explore some ideas for 
recovering structural information of the progenitors us- 
ing only a few observed HENs. Such recovery may not 
be conclusive for every detected source, but for specific 
detections where the energy and timing of detected neu- 
trinos are favorable, it can provide important information 
about the source. 



1. HEN timing prior to time of jet breakout 

The time difference between the onset of the (energy- 
dependent) observable HEN emission and the jet break- 
out can be used to constrain the possible progenitors, 
potentially even with one detected HEN. If the observa- 
tion time of a detected neutrino and the time of the jet 
breakout differs more than what is allowed for a progen- 
itor model, the respective progenitor model can be ruled 
out or weakened. For example, if one detects a HEN with 
reconstructed energy e v ss 10 3 GeV approximately ~5 sec 
before the jet breakout, Figure 4 (right) suggests that 
one can practically rule out models {I5S,I6T,40S,75S}, 
while models {15L,I5Lc} are possible progenitors. We 
note that it may be difficult to directly observe the time 
of jet breakout as the gamma-ray photosphere typically 
lies above the stellar surface. 



2. HEN relative timing 

It is also possible to constrain progenitors even if the 
jet does not break out of the star. In this case, the time 
difference between at least two observed neutrinos can 



be used to constrain possible progenitors by determin- 
ing whether the observed time difference is possible or 
likely for a given progenitor model and the reconstructed 
neutrino energies. For example, the observation of two 
neutrinos, both with energies e„ ps 10 3 GeV with ~10 sec 
time difference would rule out all models but 15L out of 
those considered in Figure 4, indicating the presence of a 
hydrogen envelope. 



3. Jet duration vs. onset of neutrino emission 

Neutrinos above e„ > I0 3 GeV for all models but 15L 
are emitted only in a small fraction of the time the jet 
spends within the progenitor. Since it is unlikely that jets 
are fine-tuned such that they stop just before breaking 
out of the star, the detection of HENs with e v > 1 3 GeV 
from confirmed astrophysical sources with no EM coun- 
terparts would make it highly probable that the progeni- 
tors have kept their hydrogen envelope prior to explosion. 



4- Neutrinos and gravitational waves 

The coincident observation of gravitational-waves and 
HENs from a common source has far-reaching astrophys- 
ical implications [27, 90]. The relation between the ob- 
served times of arrival of gravitational waves and HENs 
from a stellar core collapse could provide information on 
the stellar structure below the shock region of > 10 10 cm 
(complementary to information from HENs and EM ra- 
diation, which map the structure at and above the shock 
region) . 

The time relativistic jets take to cross the stellar en- 
velope is likely comparable to or less than the observed 
duration of prompt gamma-ray emission (as the dura- 
tion of the outflow fed by the central engine is unlikely 
to be fine-tuned to the envelope crossing time). Given 
the observed long-GRB durations of « (10 — 100) s and 
expected stellar progenitor densities at small radii, the 
jet is only likely to be able to cross the envelope within 
« (10 — 100) s if the stellar density significantly decreases 
along the rotational axis [89]. Gravitational- wave emis- 
sion from the stellar progenitor is likely connected to the 
onset of jet propagation (core-collapse and the formation 
of an accretion disk are both potential sources of gravi- 
tational waves [91]). Consequently, comparing the time 
of arrival of gravitational wave signals and HENs from a 
GRB or supernova progenitor can provide information on 
stellar densities below the shock region, and may provide 
information on the development of jets (see [89]). 



5. Extremely high energy neutrino triggers 

Extremely high energy neutrinos of energies e„ > 
100 TeV may be emitted once the relativistic jet becomes 
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sparse enough such that neutrino production is not sup- 
pressed by strong magnetic fields or high synchrotron- 
photon density [46, 47]. The emission of such extremely 
high energy neutrinos may commence only once the jet 
has advanced substantially, likely outside the helium en- 
velope [46]. These neutrinos will therefore be produced 
only after most lower energy HENs are created while the 
jet is at lower radii (we note that the central engine of 
GRBs may show activity on a longer time scale than 
the prompt gamma-ray emission [92], therefore weaker, 
longer duration neutrino emission is also plausible). Pro- 
cesses other than internal shocks, such as the jet's inter- 
action with the interstellar medium, can further produce 
so-called ultra high energy neutrinos with energies up to 
10 10 GeV that may be of interest here [38, 93, 94]. Hence, 
it may be interesting to search for > 100 TeV neutrinos in 
coincidence with < 100 TeV neutrinos that arrived prior 
to an extremely high energy neutrino up to tens of sec- 
onds. 

While HEN detectors mainly use Earth as a shield 
from atmospheric muons, extremely high energy neutri- 
nos with e„ > 1 PeV are also absorbed before they can 
cross Earth, and thus can only be detected from downgo- 
ing and horizontal directions [73]. Even though there is 
an abundant atmospheric muon background from these 
directions in the detector, these muons have lower ener- 
gies; hence, astrophysical extremely high energy neutri- 
nos arc relatively easy to identify. 

Based on the above model, a search could be performed 
for extremely high energy neutrinos in coincidence with 
< 100 TeV neutrinos that arrived from the same direc- 
tion, prior to the extremely high energy neutrino up to 
tens of seconds. This search could be particularly in- 
teresting, as < 100 TeV neutrino data have not been 
utilized in searches for astrophysical neutrinos (never- 
theless, these downgoing events are recorded and stored 
for IceCube [95]). 



VI. CONCLUSIONS 

We investigated the opacity of massive pre-supernova 
stars to high energy neutrinos (HENs) in order to ad- 
dress the question: What do the times of arrival of de- 
tected high energy neutrinos tell us about the properties 
of their source? We investigated the effect of opacity 
on the observable HEN emission from both successful 
and choked jets. In particular, we have examined vari- 
ous zero-mass main sequence (ZAMS) stellar masses from 
(12-75) M for low-metallicity non-rotating and stellar- 
metallicity rotating cases. 

For the considered progenitors, we presented the 
energy-dependent critical radius from which HENs can- 
not escape. We found that the presence of the stellar 
hydrogen envelope has a negligible effect on the optical 
depth for neutrino energies of e„ < 10 5 GeV, i.e. the 
most relevant energy range for HEN detection. The crit- 
ical radius, however, largely varies with e u within the 



helium core, which has relevant consequences on obser- 
vations. 

The neutrino emission spectrum changes as the rela- 
tivistic jet advances in the star. Considering mildly rel- 
ativistic jets (r^ w 10) and HEN production in internal 
shocks, the energy dependence of the onset of neutrino 
emission is shown in Figure 4. Such time dependence can 
provide important information on the stellar structure. 
For instance, with observation of multiple HENs, the de- 
tected neutrino energies and times provide constraints on 
stellar density at different depths in the star. The energy 
of a detected precursor HEN can also provide informa- 
tion on the maximum time frame in which the relativistic 
jet will break out of the star and become observable. 

We examined how neutrino interaction with dense stel- 
lar matter can be used to probe stellar progenitors. We 
investigated both the strong and weak-signal limits, i.e. 
when many or only a few neutrinos are detected, respec- 
tively We demonstrated that under favorable conditions 
one can use the time difference between a precursor neu- 
trino and the jet breakout to exclude some progenitor 
models. The relative times of arrival for multiple neutri- 
nos may also be sufficient to exclude progenitor models, 
even with no observed electromagnetic emission. Ad- 
ditionally, the detection of HENs with energies above 
> 10 3 GeV from choked GRBs makes it likely that the 
progenitor possessed a hydrogen envelope prior to explo- 
sion. Also, while the detection of HENs and electromag- 
netic signals can provide information on the stellar region 
at and above the shock radius, the coincident detection of 
HENs with gravitational waves [90] may be informative 
w.r.t. jet development and propagation below the shock 
radius. We proposed the use of extremely- high energy 
neutrinos (e u > 100 TeV) detected by cubic-kilometer 
neutrino detectors such as IceCube for searches for coinci- 
dent downgoing neutrinos. These extremely-high energy 
events probably arrive after other lower energy events 
that can be used in a coincident analysis. Downgoing 
neutrinos with lower energies are typically not used in 
searches due to the large atmospheric muon background 
from these directions and energies. 

A future extension of this work will be the calculation 
of neutrino fluxes from different radii, similar to the cal- 
culations of Razzaque et al. [47] , who estimated the flux 
for two different radii, but without timing information. 
Such addition to the temporal structure of neutrino en- 
ergies can provide a more detailed picture of not only the 
information in neutrinos about the progenitor, but also 
the likelihood of detecting neutrinos with given informa- 
tion content. 

We note here that the uncertainty in the reconstructed 
energy and timing of HENs introduces uncertainty in 
the measurement of the onset of HEN emissions at the 
energies of the detected neutrinos. These uncertainties 
need to be taken into account when comparing emis- 
sion models to observations. Additionally, there could 
be other factors, e.g., physics related to jet propagation, 
that we have treated schematically in this work, but could 
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have similar impact on the temporal structure of HEN 
events. However, we have pointed out some generic fea- 
tures which should motivate future work that investigates 
experimental detectability of these features. 

The results presented above aim at describing the in- 
terpretation of a set of detected high energy neutrinos 
from a collapsar event. One of the advantages of such in- 
terpretation is that it can be done independently of the 
emitted high energy neutrino spectrum inside the source, 
which greatly simplifies the understanding of observa- 
tions. An important direction of extending this work will 
be the calculation of the emitted neutrino spectrum in- 
side the source. Obtaining this time-dependent spectrum 
would let us understand the probabilities of obtaining a 
set of neutrinos with specific energies and arrival times, 
which can be used to further refine the differentiation be- 
tween progenitor scenarios. The calculation of emission 
spectra can be done similarly to the work of Razzaque et 
al. [47] , who calculated the observable spectrum for two 
specific radii for one stellar model. For accurate flux esti- 



mates, one must also take into account neutrino flavor os- 
cillations in the stellar envelope, which are non-negligible 
for the relevant energy range [67]. For precisely assessing 
the capability of differentiating between various stellar 
models, one will further need to compare the estimated 
source flux with the atmospheric neutrino background of 
neutrino detectors (see, e.g., [79, 96, 97]). 
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